Survival and residence times in disordered chains with bias 
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We present a unified framework for first-passage time and residence time of random walks in 
finite one- dimensional disordered biased systems. The derivation is based on the exact expansion of 
the backward master equation in cumulants. The dependence on the initial condition, system size, 
and bias strength is explicitly studied for models with weak and strong disorders. Application to 
thermally activated processes is also developed. 
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I. INTRODUCTION 

A large number of physical properties of diffusion and 
hopping transport of classical particles (or excitations) 
in disordered media have been investigated by means of 
random walks (RW) in disordered lattices Q ||, ||. Of 
particular interest are the effects of the finite size and 
boundary conditions on the domain of diffusion. The ab- 
sorbing boundary approach allows us to analyze when 
a process first reaches a given threshold value Q. This 
question arises in many situations and is equivalent to 
regarding the relation between the underlying dynam- 
ics of randomly evolving systems and the statisticsof ex- 
treme events for such systems. These statistics are im- 
portant in a variety of problems in engineering and ap- 
plied physics pj . Extreme phenomena are experimentally 
accessible and enable us to know the parameters of the 
stochastic dynamics. 

One quantity that naturally arises in this context is the 
time for which the particle survives before its absorption 
in the boundary sinks, i.e., the first-passage time (FPT). 
This time depends on the realization of the RW, thus 
being a random variable. The mean first-passage time 
(MFPT) is of fundamental importance for diffusion in- 
fluenced reactions, as it measures the (reciprocal) reac- 
tion rate constant. First-passage problems appear in a 
wide range of applications Jf and have a long history || . 
Recently, the fact that the MFPT is exactly equal to the 
inverse of the associated Kramers escape rate was proved 
for arbitrary time-homogeneous stochastic processes . 
Finite size effects also appear when we consider unre- 
stricted diffusion (no boundary, or boundary condition 
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too far away), however we ask for the time spent by 
the diffusing particle in finite domains. This quantity 
is the random variable known as residence time (RT). 
Unlike the FPT, which reckons the lifetime of particles 
that never abandoned a given domain, the RT involves 
the case when the particle can exit from and enter into 
the domain an unrestricted number of times. We must 
stress that unfortunately the words residence and sur- 
vival have sometimes been used as synonyms. We dis- 
tinguish these terms from the fact that the particle can 
return or not (absorption) to the interval of interest. The 
mean residence time (MRT) has importance for diffusion 
influenced catalytic reactions where reactants are local- 
ized in a finite domain of the catalyzer diffusion region. 
Experimental techniques generally called single-molecule 
spectroscopy, allow one to follow the evolution in time of 
the state of a single molecule that undergoes a confor- 
mational change (isomerization reaction). This fact has 
recently been addressed by the MRT study of a single 
sojourn in each of the states of the molecule M . The im- 
portant feature of this class of reactions is that these can 
be regarded as being one dimensional. The MRT has 
found several other applications || |l0|]. However, the 
general problem of RT distribution in random media has 
been little discussed in the physics literature 111]]. One 
of our goals is to present the FPT and RT statistics in 
random media in a unified formalism. 

The effect of bias on FPT in disordered systems has 
received attention since the first works on the survival 
fraction of particles in media with randomly distributed 
perfect tramps Jl^| . When the bias is switched on 
(through external fields), the system undergoes a sub- 
stantial change in its dynamics however small the field is. 
Several results have been reported on survival probabil- 
ity (or related quantities) for one-dimensional RW with 
disorder and bias. The models of disorder often regarded 
are the random traps and Sinai's model. Trapping in one 
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dimension is a model of strong disorder that allows exact 
results with large bias ||l3| . Sinai's model jljj is a time 
discrete RW in one dimension with asymmetrical transi- 
tion probabilities that fulfill a certain condition. In this 
way, in the Sinai's model the disorder is coupled with 
the bias strength. The FPT problem for Sinai's model 
has been extensively studied |15fl. In this work we con- 
sider a RW on a chain with site disorder in the presence 
of global (site independent) bias. We analyze weak and 
strong situations of disorder [[[} |l6| . For weak disorder all 
the inverse moments of the distribution of RW hopping 
rates are finite, whereas all these moments diverge if the 
disorder is strong. 

A successful theory for FPT statistics in disordered 
media is the finite effective medium approximation 
(FEMA) FEMA combines an exact expansion of 

the survival probability equation in a disordered medium, 
with the effective medium approximation ]l7j] . This 
scheme allows a perturbative analysis around the effec- 
tive homogeneous medium in the long time limit for 
weak and strong disordered models. FEMA also pro- 
vides a self-consistent truncation criterion. The exten- 
sion of FEMA to biased media was presented in Ref . [jl8| , 
where we got perturbative expansions for small bias and 
weak disorder. In Ref. fll9|] , another extension of FEMA 
was carried out for periodically forced boundary condi- 
tions. In the present paper, we obtain, in the guidelines 
of the expansions developed in FEMA, the exact equa- 
tions for MFPT and MRT and construct their solutions 
in the leading relevant order for small bias. MRT for 
one-dimensional diffusion in a constant field and biased 
chains was analyzed in Ref. However, we could not 
find explicit expressions for the RT distribution in disor- 
dered media in the literature. Therefore, another goal of 
our work is to consider the mixed effects of disorder and 
bias in the FPT and RT distributions. 



This paper presents the survival and residence times 
statistics. In Sec. II we define the survival and resi- 
dence probabilities and construct their expressions from 
the conditional probability of the random walk on a 
chain. The random biased model is described in Sec. p| , 
whereas in Sec. [V the homogeneous (nondisordered) 



chain with bias is treated analytically. In Sec. [v| we in- 
troduce the projection operator to average over disorder, 
obtaining the main equations. The weak disordered case 
is analyzed in Sec. VI and strong disor dered cases are 
considered in Sec. VII. Finally, in Sec. VI II thermally 



activated processes are considered and Sec. LX provides 
the concluding remarks. The mathematical details of the 
paper are condensed in two appendixes. In Appendix |a] 
the survival and residence probabilities for homogeneous 
chains are exactly calculated, in Appendix [b| we study 
the Green's functions in the presence of bias, and in 
Appendix P we evaluate the relevant cumulants used in 



II. SURVIVAL AND RESIDENCE TIMES 
STATISTICS 

The dynamical behavior of random one-dimensional 
systems can be described by the one-step master equa- 
tion § 

d t P(n,t\n ,t ) = P(n - l,t\n ,t Q ) 

+ w- +1 P{n + l,t\n ,t ) (2.1) 
- (w+ + w n )P(n,t\n ,t ) , 

where Wn ' is the transition probability per unit time 
from site n to n+1 (n— 1 ). P(n, t\no, to) is the conditional 
probability of finding the walker at site n at time i, given 
that it was at site no at time to (< t) and a particular 
configuration of {w^;}. We assume that {u>^} is a set of 
positive independent identically distributed random vari- 
ables. In this description, the disorder is modeled by the 
distribution, p(w), assigned to these random variables. 
For a given realization of {«^} (quenched disorder) we 
get a M arkovian stochastic dynamics. We can write 
Eq. (2.1) in matrix notation: d t P(t\to) = H P(t|i ), 
where 



n+1 n' 



(2.2) 



and P(tol^o) = I- Thus, the formal solution of Eq. (2.1) 
is P(£|£o) = exp((i — io)H). This solution obeys the 
backward master equation too, <9 t P(t|to) = ~P(t\to) H, 
for the same initial condition P0[ . 

In this work we consider a RW on a chain and we 
address the question about the survival and residence 
probabilities in the finite interval D = [— L, L\. The first 
is the probability, S no (t\to), of remaining in D (without 
exiting) at time t if the walker initially began at site no € 
D. The second, instead, is the probability, R no (t\to), of 
finding the particle within the domain D at time t, given 
that it initially began at site no (not necessary in D). 
Therefore, the residence probability is defined by 



Rn (t\t Q ) = ^2 P(n,t\n ,t ) , 



(2.3) 



n£D 



where P(n, t\no, to) is the solution of Eq. (2.1). Due to 
the fact that the master equation links the probabilities 
for all sites of the chain, the residence problem involves 
an infinite matrix. 

To compute the survival probability, we need to elim- 
inate contributions from trajectories returning to the in- 
terval D after having left it. To do it we must find the so- 



lution Po{n,t\no,to) of Eq. (2^) with absorbing bound- 
aries in the interval's extremes |4|. Thus, Po{n, t\no, to) 
results in the solution of dtPD(t\to) = Hd Po(t\to), 
where 
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H, 



if n' e D 
otherwise 



(2.4) 
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Hence, the survival probability results in 



Sn (*|*o) 



neD 



P D (n,t\n ,t ) . 



(2.5) 



This definition and the fact that for the survival problem, 
we only need to consider no € D allow us to work with a 
finite square matrix of dimension N x N, N = 2L + 1 
being the number of sites in D. 

The presented view of the survival problem is adequate 
for chains with a fixed number of sites. Nevertheless, we 
are interested in general expressions for domains with an 
arbitrary number TV of sites and we want to work out 
survival and residence problems simultaneously. Due to 
the time-homogeneous invariance of the problems we take 
t = from here and throughout the rest of the work. 
We now consider the vector function F(t) whose com- 
ponents are F no (t) = Yln&D P{n-,t\no, 0). The evolution 
equation for this function follows from the backward mas- 
ter equation and results in dtF(t) = F( t), w here 
is the transpose matrix of H. Using Eq. ( [2.2| ), we can 
write the last equation in components 



d t F 



n (t) = 
+ 



< [Fn + l(t) - F no (t)} 

w~ [F no -i{t) - F no (t)} 



(2.6) 



Thus, the residence probability is the solution of Eq. (|2 
with the initial condition, 



Rn (t -0)-\ q oth ° rwise 



(2.7) 



and boundary condition at infinity: R no (t) — > for 
| no | — > oo for all finite t. On the other hand, the sur- 
vival pro babi lity is the solution of the generic adjoint 
equation ( [2.6] ) for the infinite chain with the initial con- 
dition Sn (t = 0) — 1 for all no S D. Here, the artificial 
boundary condition S no (t) = for all t if no = — (L + 1) 
or no = L + 1 must be used to prevent the back flow of 
the probability into the interval |2lJ . 

The survival probability decreases monotonically in 
time from unity to zero. Let us now introduce the first- 
passage time distribution (FPTD) f no (t), i.e., the prob- 
ability density of exit D at a time between t and t + dt; 
then f no (t) = -d t S„ (t) @. The MFPT is the first mo- 
ment T no (if it exists) of FPTD, 



tf n „(t)dt 



If t S no (t) — > for t — > oo , then 



Tn — / S ng (i) dt 
Jo 



(2- 



(2.9) 



The residence probability does not necessarily decrease to 
zero at infinitely long times. Moreover, it need not even 
be monotonic in time. Thus, the residence time density 
is generally not equal to the negative time derivative of 



the residence probability |Tl|]. Nevertheless, we can de- 
fine the MKT, r„ , in a manner analogous to Eq. ([2.9]), 
namely, 



Rn (t)dt 



(2.10) 



Thus, from Eqs. (gj) and Q2.10[ ), we obtain MFPT and 
MRT from the asymptotic limit of the Laplace trans- 
formed (denoted by hats) survival and residence proba- 
bilities |E2|, 



lim S no (s) , 



limits) . 

s— >0 



(2.11a) 
(2.11b) 



For MFPT, this limit exists if S m (s) ~ T no +cs x , where 
c and T no are assumed constants and x > 0. In this 
manner, from f no (s) = 1 — sS no (s), the normalization 
condition of the FPTD is also guaranteed: J °° /„„ (t) dt = 

lim^o fn (s) = 1. 



III. RANDOM BIASED MODEL 

In this work, we are interested in the interplay between 
the bias and the disorder in the transition probabilities. 
For this goal we take 



in 



(3.1) 



where a and b are positive constants and are taken 
to be independent but identically distributed random 
variables with (£„) = 0. This form of jump transitions 
involves an ordered biased background with a superim- 
posed random medium. The strength of the bias is given 
by the ratio between a and b and the disorder is char- 
acterized by the distribution of variables Without 
lost of generality we assume a > b and in consequence we 
must to impose the restriction > — b. This restriction 
guarantees the positivity of jump probabilities {w^}. We 
introduce the parameter e for bias strength by b/a = 1 — e 
and we take < e < 1, so that the bias field points to 
the right. This election of parameters allows us to fo- 
cus our attention in the small bias limit and to st udy 
the transition to the symmetric diffusive behavior p3| . 
The Laplace transform of th e evolution equation for our 
model results from Eq. ([2.6D, 



sF n (s) - F n (t = 0) = [K. b +CnX°] F n (s) , (3.2) 
where we have introduced the operators 



K° = £++£- -21 , 

JC b = a (S+ -l) + b (S- -T) . 



(3.3) 



E^- are shifting operators (£^g n = 9n±i) and X is the 
identity operator. Equation (|3.2|) must be solved with 
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the boundary conditions corresponding to each problem, 



S-(l+i)(s) = S (L+1) (s) = 0, 
Rn(s) — > for \tl\ — > oo , 

V s > 0. The initial condition is given by 
F n (t = 0) = l Q otherwise _ 



(3.4a) 
(3.4b) 



(3.5) 



Remember that for FPT we only need to consider n E D. 

The classes of disorder analyzed are generalizations of 
standard cases in the literature jl], Our expressions 
are constructed introducing the parameter e in such a 
way that we guarantee the positivity of transition rates 
and reproduce the known expressions in the limit e — > 0. 
We have considered the following three classes of disorder 
for the transfer rate w = wt ■ 



The expressions given in Eqs. ( p.6[ ) and (3.7) can be 
obtained from the corresponding expressions given in 
Eqs. and in the limit a -> 0. Class (g) cor- 

responds to the situation of weak disorder. There, the 
mean-square displacement of the RW behaves like t for 
long times. Classes ([b]) and (£]) become strong disordered 
cases for e — > 0, and correspond to situations of anoma- 
lous diffusion. For class (|]), fix oc ln(l/e) and /3* oc e 1_fe 
if fc > 1. For class ([J, /3fc oc e 1 ^ k ^ a . In the strong dis- 
order limit, (n(t) 2 ) behaves for long times as t/hit and 
t 2(i-a)/(2-o) ) f or cases Q and (|) ) respectively. Though, 
the MFPT in the presence of strong disorder is a diver- 
gent quantity p6[ . We will show in Sec. VII that our 
model allows us to study the transition to strong disor- 
der in the limit e going to zero, i.e., the zero bias limit. 



IV. NONDISORDERED BIASED CHAIN 



a. The mean values of the inverse moments of jump 
transition w, (3k = ((l/w) fc ), are finite quantities 
for all k > 1, e > and remain finite in the limit 
e -► 0. 

b. The probability distribution p(w) is 



(w) = { B if W 6 ( ae '^ s ) 
^ ' ] otherwise , 



(3.6) 



where the values of B and £Ib are fixed by the 
normalization condition (J Q p(w) dw = 1) and the 
fact that (w) = a, 



B = (2(l-e)a) 1 
tt B = (2- e)a . 



(3.7a) 
(3.7b) 



c. The probability distribution p(w) is 
p(w) ■■ 



Cw- 




if w G (d£, Oc) 
otherwise , 



(3.8) 



where < a < 1 and the values of C and Sic* 
are also fixed by the normalization condition and 
(w) = a. For small e it gives 



C 



(I -a) 



2-a 



[{2 -a) a] 



l + (2-a) Q ((l-^e) 1 



(3.9a) 



The nondis ordered case is obtained from the trivial 
distribution p(w) — 5(w — a). Therefore, basic results 
about the survival and residence probabilities in nondis- 
ordered chains can be easily obtained from the equation: 
sF n (s) — F„(0) = JC b F n (s), with the boundary condi- 
tions given by Eq. ([O]). In particular, exact expressions 
for MFPT and MRT for a homogeneous biased chain are 
given by (see Appendix ^ for detailed calculations) 



L + l-n 2 (L + 1) 7 r ' 



a(l-7) a(l - 7) 7 -( i + 1 ) - 7 L +! 



(4.1) 



with — L < n < L, and 



r 2L+ 1 



1 

- < 

a 



1-7 

L-n 1 - 7 ™+ i + 1 

1^7 + a-7) 2 

1 - 7 2i+1 
[ (1-7) 2 7 



n < -L 

-L <n< L (4.2) 
n > L , 



where 7 = 6/a < 1. Figure [l] shows the behavior of r n 
for some values of 7. We would like to stress that for 
n < —L, MRT is a constant proportional to the width of 
the interval (2 L + l), whereas for n — > 00, MRT vanishes. 
Given that the bias points to the right, for any initial con- 
dition at the left of the interval of interest, MRT is equal 
to the transit time across the interval. For a given initial 
condition, MRT is greater whereas th e bi as is smaller. 
For one way motion (7 = 0), from Eq. (4.2) we obtain 



c 



a 



(3.9b) 




(4.3) 
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FIG. 1: MRT for a nondisordered chain [as given by Eq. ([4.2| )1 
plotted against the discrete initial condition n, with a — 1 
and L — 10. The solid lines are only to guide the eye. 



whereas in the small bias limit, i.e., 7 = 1 — e, e << 1, 
taking a constant results in 



2L + 1 



a e 



1 

_ (L + n) (L + n + 1) 



1 — n e 



2(2Z + 1) 



, n < — L 
, ne D 
, n > L . 



(4.4) 



It is worthwhile to emphasize that MFPT exhibits 
a crossover from the drift (strong bias) regime [T n = 
(L + 1 — n-)/a)] to the diffusive (small bias) regime 
(T n oc ((L + l) 2 — n 2 )/2a). The diffusive behavior is also 
present for MRT for n S D, but it is not in the leading 
term. Moreover, MFPT remains finite in the limit e — > 
for finite dom ains, whereas MRT diverges as we can see 
from Eq. (14). Thus, the MRT is not a defined quantity 
for unbiased chains pi| . Expressions for the MFPT in 
nondisordered biased chains following from Eq. ( |4.l| ) , and 
the study of the drift and diffusive regimes were reported 
inRef.Pl. 



V. PROJECTION OPERATOR AVERAGE 

The basic equation f or t he evolution of probabilities 
has been written in Eq. (|]^) . The operators in this equa- 
tion explicitly show the splitting of the transition prob- 
ability in an average biased part ((w„) — a, (w~) = b) 
and a random nonbiased part ( £„). Defining the operator 
A = KP, let us rewrite Eq. (|3.2|) as 



sF(s)-F(0)= [K, b + A] F(s) 



(5.1) 



Our goal in this section is to obtain exact equations for 
the averaged survival and residence probabilities. This 
average can be formally carried out introducing a pro- 
jection operator V (V 2 = V) that averages over the 
joint probability density of variables {£, n }- (F) =VF, 



F = (F) + (l-T) F. Applying the operator V to Eq. fl5.l]) 
we obtain 

s (F) - F(0) = K. b {F) + VA(F) + VA(l - V)F . (5.2) 



Also, applying the operator (1 — V) to Eq. ( p.l[ ) we arrive 
to 



s{l-V)F = JC h {\ -V)F+ (1 -V)A(F) 
+ (1 - *P)A(1 - V)F . 



(5.3) 



A formal solution of Eq. ( p.3[ ) can be obtained using the 
Green's function for the nondisordered chain, 



goo 



(5.4) 



Appl ying G to Eq. ( p.3| ) and using the definition given in 
Eq. (5.4), results in 



(1 - V)F = G (1 - T)A(F) + (1 - V)A(1 - V)F 



(5.5) 



Equation (p.5|) can be iteratively solved for (1 — V)F, 



OC t 

(i-v)f = J2[g(i-v)a\ (F) 



k=l 



(5.6) 



Putting this formal solution in Eq. (5^) we find a closed 
exact equation for the average probability (F) , 



s(F)-F(0)=lC b (F) + (j2[AG(l-V) 

\k=0 



A J (F) . 
(5.7) 

The operator G is solution of the equation, 

(s-IC b )G = t, (5.8) 

with the boundary conditions 

G£( L+ i )m (s) = Gf L+1)m (s) = Vs and m eD , 

(5.9a) 

G„ m (s) — * for \n\ — * oo and m finite , (5.9b) 

where the superscript S (R) corresponds to the survival 
(residence) problem. Exact expressions of the Green's 
functions are calculated in the Appendix M 

We will find it useful to write Eq. ( |5?7|) in compo- 
nents. For this task, we use the explicit form of A, 
the Terwiel's cumulants [^5| of the random variables 
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and we define the propagator: J nm (s) = fC° G nm (s), 
where the operator K? acts on the first index of G nm {s). 



Expl icit e xpressions for J nm (s) are given by Eqs. flBTj ) 
and (B12) in Appendix H. In this manner we can write 



s(F n )-F n (0) = JC b (F n )+Y, E 

p=0 mi,... ,m p 
X (£n£mi • • ■ £,m p )T Jnmi (s) 
X Jm 1 m 2 ( s ) • ■ ■ ^m p _im p (s) /C (F mp ) . 

(5.10) 

We must understand that rn = n : an< i f° r the survival 
problem we have the additional restriction mi, . . . , to p S 
Z). The exact effective backward equation given by 
Eq. (5.10) can be rewritten as 



s(F n )-F n (0) = JC\F n 



E E 



p=0 



m p ^m p _ 1 

X (ipnlpmi ■ ■ ■ V'mjT Jnmi (s) 

X imims ( 5 ) • • • Jm p -\rn p { s ) K-> (F mp ) , 

(5.11) 

where we have used the definition of Terwiel's cumulants. 
Here, we have summed up all the terms containing the 
diagonal parts of J nm (s) through the introduction of the 
random operator tpk(s) defined by 



ipk(s) 



OO 

E 

ik=0 



[Jkk(s)i k (l-V)} lk & 



(5.12) 



This operator acts on any disorder-dependent quantity 
at its right. The geometrical sum in Eq. (5.12) can be 
evaluated, resulting in 

Ms) = M fc (*) - - M t S V^\ \ VM k ( S ) , (5.13) 



I + (Mk(s)) Jkk(s) 



where 



M k (s) 



1 — Cfe Jkk(s) 



(5.14) 



In the following, we ta ke the limit s — ► in order 
to obtain from Eq. ( 5.11 ) the corresponding equations 
for the MFPT an d MK T in disordered media, which are 
defined from Eq. ( 2.11 ) by 



(T„) = lim (S n (s) 

S— HJ \ 

(r n ) = lim {Rn(s) 

s— >0 \ 



(5.15a) 
(5.15b) 



In this limit the propagator J nm (s) can be written as 
A nm + O(s), and the exact expression for A nm in the 
FPT (RT) problem is given by Eq. flBl^ ) [Eq. ( p!5| )] in 



Appendix |b[ Therefore, the resulting equation for the 
averaged MFTP ((T n )) is 



S n (t = 0) = JC\T n ) + E E ■ ■ ■ *!>■ 



ni. p /T 



p=0 m,^n 



X A-3 A S 

nmi mi m2 



A 5 K° IT ) 

■ ll m p - 1 m p /v - \- L m p l 



(5.16) 

Its solution must satisfy the boundary conditions: 
(T_(£_|_i)) = = 0. A similar equ ation for the 

averaged MRT ((r„)) is obtained from Eq. (5.16) replac- 
ing the S quantities by the corresponding R quantities 
and imposing the boundary conditions: (r„) = const for 
n < ~L and (t„) — > for n — ► +oo (given that the bias 
field points to right). 

We, additionally, consider the case of small bias. Thus, 
the expressions for A„ m can be further expanded, and 
taking a constant results in 

1 ' -(L 



A 



2(L + l)a L 



2{L + \)a 



I 2{L + l)a 



~(L- 



l)]e 



A, 



,71 < TO 

— 1 , n = to 
e , ft > to 



, 71 < 771 
, 71 = 771 
, 71 > 771 , 

(5.17a) 
(5.17b) 



Again, the superscript S (R) corresponds to the sur- 
vival (residence) pro blem. Hence, A„ m oc e for n ^ m 
and, from Eq. (B15), the diagonal components of the 
propagator in the residence problem result independent 
of e. 



VI. WEAK DISORDER 



For disorder of class (0), the quantities 0k are fi- 



nite and we obtain that (VViVv, 



oc e 



and 



(V'nV'mi • ■ - ipm^T is independent of e. Thus, Eq. ( 5.16| ) 
is a perturbative expansion in the sense that (T n ) can be 
calculated up to order e q truncating the p series in the q 
term (p = 0, . . . , q). The corresponding equation for (r„) 
is strictly a perturbative expansion given that the con- 
tribution to the order e q comes entirely from the term 
p = q. 

In Eq. ( [3.16 ) , the cumulant for p — 1 consists of two 
independent random operators, so it vanishes (see Ap- 
pendix ^J). Therefore, it turns out that only the t erm 
with p — contributes to order e. From Eq. (C14) re- 
sults 



WOt = l 3 i 1 -a-aT 



1—71 



2(L+1) 



(6.1) 
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where we have introduced the fluctuation of the quenched 
disorder: T — {fii — (3i)/(3 2 . Up to order e, taking a 
constant and usin g JC b ss a [lC° — (£~ — I) e] , the explicit 
form of Eq. ( |5.16|) is 



1 = (3^IC (T n ) 
L + l 



T 



2 (L+l) 



K? + {£- - 1) 



(Tn) 



0(e 2 ). 
(6.2) 



For e = 0, Eq. ( |6 . 2| ) immediately gives the well- 
known MFPT for the unbiased case: (T n (e = 0)) = 
[(L + l) 2 — n 2 ] / (2 /3f 1 ), where we can see that the ef- 
fect of weak disorder is to replace the constant a by the 
effective coefficient (5^ 1 . To construct the consistent solu- 
tion up to order e of Eq. (|6.2]), we propose the expression 



(Tn) = 



(L + l) 



2(3- 



[l + {An + B)e] 



(6.3) 



which immediately satisfies the boundary conditions: 
(T-( L +i)) = (Tl+i) = 0. To fit the constants A and 



B we substitute this expression in Eq. (S.2) and retain 
only the t erm s up to order e. If the factor of e in the ex- 
pression ( |6.3|) were a polynomial of degree greater than 
1, we can easily see that the coefficients of the terms of 
degree greater than 1 vanish. Thus, we obtain for n € D 



(T n 



(L + lf 



3-2n L + l-n/3 
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2 (L + l) 



F I a(3\ e 



(6.4) 



From this expression we can analyze the interplay be- 
tween the bias and the fluctuation of the quenched disor- 
der T . In the presence of bias, the escape time from the 
finite interval increases for all initial conditions n, with 
respect to the unbiased case, if the fluctuation is large 
enough. This is an important result related to the con- 
trol of the trapping process. This fact was reported in 
Ref. [[l8] where the solution of Eq. ( |6 . 2| ) was constructed 
from FEMA. Now, we can evaluate the difference be- 
tween the exact solution (order e) given by Eq. (6.4) 
and that approximation, obtaining that the correction 



to FEMA is 



(L + l) 2 



6 (L+l) 



nJ-ae. 



In the residence problem, from Eqs. (C9) and (C14) re- 
sults (VVi)t = Pi 1 ~ a - Thus, the corresponding equation 
for the averaged MRT is 

-R n (0) = (3^ K° (r„) — ae (£~ — I) (r„) + 0(e 2 ) , 

(6.5) 

where i? n (0) is given by Eq. ( |3.5| ). Now we have to deal 
with a backward master equation with constant coeffi- 
cients. In this case, the fluctuations of disorder are not 
present. Formally, this equation is equal to the one cor- 
responding to a nondisordered chain in the small limit 



bias. Therefore, from Eq. (4.4), with the substitutions 



a ft" 1 , 

b ^(l-ojSie), 
e — > a Pie, 

we obtain, up to order e, the solution 

21 - ! I (L + n)(L + n + l) 
T n w < 1 - T^rh -afoc ,n € D 



(6.6) 



, n < —L 



2 (2 L+l) 

1 — n a (3\ e 



VII. STRONG DISORDER 



,n > L . 
(6.7) 



For the classes of strong disorder, (3k are divergent 
quantities in the limit e — > 0, 



A«|lne|/(2a) 
P^M^ (*>1) 



2a(k - 1) 



, , for class (|b|), 



C_ 

■ 1 + a 



(ae) 



1—k—oi 



for class 



(7.1a) 



(7.1b) 



where < a < 1 and C is given in Eq. (|]^) . Also in this 
limit, Terwiel's cumulants diverge except the first one, 
(i>n)T, which vanishes. However, all the terms of the p 
series in the Eq. fl5.16| ) vanish and the one corresponding 
to p = is the leading term. We reckon the relevant 
cumulants in Appendix 

For a disorder of class (|]) we obtain T ~ A , and 



from Eq.(6.1) results (V>n)f ~ (3\ X ■ Therefore, the lead- 
ing terms in the equation for the averaged MFPT result 
in 



-\=^ X K? (T n )+0 



In el 



(7.2) 



This is a backward master equation for the unbiased 
chain with constant coefficients whose solution is given 
by (T n ) = [(L + l) 2 -™ 2 ] /(2/37/ 1 ). Hence, the aver- 
aged MFPT diverges as (3\ oc | In ej. 

On the other hand, for a disorder of class (Q) we obtain 
T w ^(ae)"- 1 , and from Eq. (|+]) results (^„)f w 

Pi 1 {} ~ l+S 2(£+i) ) ' Therefore, the MFPT's equation 
is 



-1=1- 



L + l 



1 + a 2 (L+l) 



(3^K. U (T n )+0{e) 



(7.3) 



In this case, we get an unbiased backward master equa- 
tion with linear coefficients. General solutions of this 
equation will be given elsewhere. Nevertheless, we hold 
our attention in the divergent behavior of the averaged 



8 




FIG. 2: Disordered potential in the chain, d is the distance 
between nearest-neighbor sites. The particle jumps to neigh- 
boring sites through potential barriers equal to 4> — (f>„, where 
4> n are independent identically distributed random variables 
with zero mean. 



0. From Eq. (7.2) we can see that 



MFPT for e 
(T n ) ocft cx e~ a . 

FEMA consists in introducing an effect ive n on- 
homogeneous medium T n and truncating Eq. ( 5.1 6| ) to 
the first term: —I = JC b (T n ). Here, the operator K, b 
is constructed with the constants a + T n and b + T n . 
The effective rate, T n , is fixed imposing the condition 
(ipn(r n ))T = 0- The resulting solutions give the ex- 
act laws: (T n ) cx |lne|, for a disorder of class (|b|) and 
(T n ) cx e~ Q , for a disorder of class ([]), however in the last 
case the predicted n-dependent coefficient is not exact. 
Strikingly, the predicted divergence laws for the averaged 
MFPT in biased chains agree with the corresponding laws 
for the survival probability, obtained for the unbiased 
case in the limit s — > Jl6| . The behavior of the survival 
probability obtained from FEMA is (S n (s)) cx | lns|, for 
a disorder of class (^), and (S n (s)) cx s _Q , for a disorder 
of class (|). Accordingly, both problems have the same 
exponents despite their different nature. 

In the residence problem, for both classes of strong 
disorder we get that (if^n)^ ~ Pi > an d then the equation 
for the averaged MRT results ~i?„(0) = /3f 1 K? (r„), i.e., 
an unbiased equation. Therefore, the averaged MRT is 
not a defined quantity for strong disorder. 



VIII. THERMALLY ACTIVATED PROCESSES 

Before concluding this work, we want to analyze the 
physical realization corresponding to thermally activated 
processes in weakly disordered chains. For this goal we 
suppose a particle that moves in a periodic potential U 
with minima separated by sharp maxima of height 0.The 
jump probability per unit time to neighbor sites involves 
the Arrhenius factor: W exp(— (3(f)) with P = (fc^T) -1 , 
where ks is Boltzmann's constant and T is the temper- 
ature. In the presence of disorder, the transition proba- 
bilities are still symmetrical but depend on the site [see 
Fig. (§)]: w n = Wexp[-f3(<p- c/) n )], where (0„) = 0, 
\<Pn\ "C 4>, and we suppose that the random potential is 



smaller than the thermal energy, i.e., P\<p n \ <C 1. Now, 
suppose that an external field E (toward right) is added, 
then the force q E is applied on the particle of "charge" q. 
The presence of E alters the heights of potential barriers. 
The new jump probabilities are nonsymmetric, 

w± = Weiq>[-(3(cj)-4> n TqEd/2)}, (8.1) 

where d is the separation between neighbor sites. In 
the low field limit, we can use the aproximation 
exp\ -p( -(f) n TqE d/2)] « 1 + f3 cp n ± q E d/2, thus [see 
Eq. @], 

a = Wcxp(-f3(j)) (1 +/3qEd/2) , (8.2a) 
b = Wexp(-/34>)(1- PqEd/2) , (8.2b) 
£„ = WeM-P0)P0n- (8.2c) 



Hence, 7 = b/a ~ 1 — j3 q E d, i.e., e = (3 q E d. 

Defining ao = W exp(— 13(f), we can write a = ao(l + 
e/2), b — a (l — e/2), and ^„ = a (3(j) n . Therefore, 
we must adapt the equations of Sec. [v| to incorporate 
the lineal dependence of a on e. In the following, it 
will be useful to define the quantities w = ao + £ n , 
Wk = (w- k ), and T = JhJJ3 ° 



1 which are inde- 



pendent of e. From Eqs. (B14) and Eqs. (B15) results 
that Eqs. (5.17) are valid with the substitution a — > a . 

Additionally, we can see that (ip n )T = Pi — «o ~ 
aoTOne, where 9 n is given by Eq. (p9|). Now, using 
that K b w a Q [JC° - (£+ - £~) e/2], the explicit form of 
the MFPT's equation up to order e is 



-l = (3 l IC°(T n )-a 



F0 n K° + ^ (£+-£-) 



e(T n ). 
(8.3) 



Proposing the solution given by Eq. ( p.3[ ) , we now obtain 
(L + lf-n 2 



(Tn) 



2 Pi 



n L + l — n/3 —\ 

3+ 2(L + 1) T ) a0/?l£ 



(8.4) 



In order to obtain the explicit dependence on tem- 
perature, we additionally must replace the disorder- 
dependent quantities according to f3 x w (1 + /3 2 (^))/ao 
and^«/3 2 (02). 

On the other hand, the MKT's equation up to order e 
results 



-i?„(0)=/3 1 1 /C° (r„) + ^e(^ 



£' 



(8.5) 



Again, this equation is equal to that corresponding to a 
nondisordered chain for small bias with transition rates 
given by the substitutions 



Ti 
Ti 1 



ao 
2 

a a 



«o Pi e , 



(8.6) 
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Thus, u p to order e, MRT is equal to the expression given 
by Eq. (3/7) with the changes a — > Qq and j3\ — > f3 1 . 



APPENDIX A: SURVIVAL AND RESIDENCE 
PROBABILITIES IN A NONDISORDERED 
BIASED CHAIN 



IX. CONCLUSIONS 

We have presented a unified framework for the FPT 
and RT statistics in finite disordered chains with bias. 
Exact equations for the quantities averaged over disor- 
der were obtained for both problems and its solutions up 
to first order in the bias parameter were constructed re- 
taining the full dependence on the system's size and the 
initial condition. 

We have studied the FPT and RT problems for three 
models of disorder. For weak disorder, the inverse mo- 
ments of the transition probabilities are finite, and we get 
that the bias becomes a control parameter for the MFPT, 
coupled with the fluctuation of the disorder. The MRT is 
only defined in the presence of bias, and for weak disorder 
the MRT's expressions are obtained from the nondisor- 
dered case renormalizing the transition constants. For 
strong disordered cases, for which the MFPT is not de- 
fined in unbiased chains, the bias allows us to study the 
divergent behavior. Amazingly, the exponents of the di- 
vergences in MFPT obtained for vanishing bias coincide 
with that obtained for the averaged survival probability 
in the long time regime. The MRT is a divergent quan- 
tity under strong disorder because the strength param- 
eter is not present in the corresponding equation in the 
small bias limit, and the MRT is not defined for unbiased 
walks. 

We complete the work with three appendixes. In Ap- 
pendix ^ the derivation of survival and residence proba- 
bilities for nondisordered chains is completely developed. 
Apppendix |b] is devoted to the detailed calculation of 
Green's function in a chain with bias. Exact expressions 
are given, which display the full dependence on the sys- 
tem's size, initial condition, and bias strength. The eval- 
uation of the relevant Terwiel's cumulants is reported in 
Appendix ^. 
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We address the question about the survival and resi- 
dence probabilities for a RW in the finite interval D = 
[-M, L] of a homogeneous chain. In the absence of dis- 
order, the Laplace transform of the evolution equation 
gives sF n (s) — F n (0) = JC b F n {s). This equation can be 
written as 

[a£+ + b£- - (s + a + b) l] F n (s) = -F n (0) , (Al) 

and must be solved with the boundary conditions corre- 
sponding to each problem, 

^_(M+i)(s)=%+i)(s) = 0, (A2a) 
R n (s) -> for \n\ -> oo , (A2b) 

Vs > 0. Defining the vector A such that A n+ i = F n (s), 
Eq. ( |Al| ) can be written as 



(a (£ 



+ \2 



(s + a + b) £+ + bl) A n = -F„(0) . (A3) 



Equation (A3) is a second-order linear difference equa- 
tion. A particular solution for the inhomogeneous equa- 
tion is the constant A n — l/s. Proposing A n = x n as 
a solution of the corresponding homogeneous equation, 
it gives x 1 — (r + 1 + 7) x + 7 = 0, where r — s/a and 
7 = b/a < 1. The roots of this second-order equation are 
given by 



1 



2:1,2 



+ 1 + 7± y/{r + 1 + 7) 2 - 4 7 



(A4) 



and satisfy x\X2 =7 with x\ > 1. Then, the general 
solution of Eq. (Al) can be written as 



F n {s) = -+c 1 x 7 { + C2x2 , ne D. (A5) 
s 

For the survival problem the constants c\ and C2 are fixed 
by imposing the boundary conditions given by Eq. ( A2a| ). 
Thus, we found for the survival probability the following 
expression: 



J 



L+M+2 



,X+M+2 



S n (s) = - I 1 - 



n-L-l 







L+M+2 



L+M+2\ 



n-L-l 



L+M+2 



L+M+2 



(A6) 



This equation is invariant under the transformation x\ <-> 
x%. From Eq. (A6) and f n (s) — 1 — sS n (s), we can 



immediately obtain the FPTD for the ordered chain. 
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For the residence problem, the boundary condi- 
tions ( A2b| ) impose that 



Rn(3) 



dx x\ for n < -M 
d 2 x 2 for n > L . 



(A7) 



Now, the const ant s c\, c 2 , d\, and d 2 are fixed by writing 
explicitly Eq. (|Al]) for the sites n = —(M + 1), -M, L, 
and L + l. Thus, we finally found the expression for the 
residence probability, 



R n {s) 



1 



s (xi - x 2 ) 



L+M+l 



1) x\- h , n < -M 



(l-x 2 ) (4 

1 - [(1 - x 2 ) x\~ L + (xi - 1) , n e D 

(xi - 1) (1 - x% +M+l ) x^ L n>L. 

(A8) 



Taking the limits of Eq. (2.11) it gives, on one hand the 
MFPT for the ordered chain, 



L + l-n L + M + 2 



7-7 



L+l 



a(l-7) a(l-7) 7~(- M + 1 ) - 7 i + 1 
with — M < n < L, and on the other hand the MRT, 

r l + m + 1 



(A9) 



1-7 

£-71 1 - 7 «+ M +! 

1^7 + (1-7) 2 

_ !. 

(l- 7 )2 7 



ra < -M 
—M <n<L 
n> L . 



Equa tions Q4.1[ ) and ( ]4.2|) follow from Eqs. (JA9J) and 
and ( A10[ ), respectively, taking M = L. 



APPENDIX B: GREEN'S FUNCTION IN A 
CHAIN WITH BIAS 

In this section we are concerned with the solution of 



written as 

aG n+ i m (s) +bG n -i m (s) - (s + a + b) G nm (s) = -8 nm . 

(Bl) 

We must solve this equation with the boundary condi- 
tions corresponding to each problem, 



G^( M +i)m( s ) = G (L+i),n( s ) — V s and m (z D , 



Gnm(s) —> for \n\ — > oo and m finite 



(B2a) 



(B2b) 



Here, the superscript S (R) denotes the solution corre- 
sponding to the surviva l (r esidence) problem. 

The solution of Eq. (Bl) in a finite domain (survival 
problem) can be constructed using the method of im- 
ages M. It consists in summing to the Green's function 
in the absence of boundaries, terms corresponding to the 
specular image of the index n with respect to the bound- 
ary considered. For absorbing boundaries the image must 
have a negative sign. Moreover, in the presence of bias, 
the image must change the bias direction (a <-» b). Addi- 
tionally, in a closed domain, each boundary reflects the 
image of the other boundary. This fact introduces an in- 
finite set of images. However, we adopt here the simpler 
algebraic approach developed in Appendix [a|. Thus, we 
propose 



(A10) d s nm (s) 



7i x'{ 
Si xl 



■ 72 x'i 
5 2 x% 



for n < m 
for n > m 



with n, m G D , 
(B3) 



where the functions x\_ 2 ai 'C defined in Eq. (A4). The 
constants 71, 72, Si, and 82 are solutions of the set of 
four algebraic equations that result from imposing the 
boundary conditions given by Eq. (B2a), imposing the 
equality between both expressions given by Eq. ( |B3| ) for 
n = m, and writing Eq. (Bl) for n = m. Additionally, 



Eq. (5.8). In components, this backward equation can be using the relation a x\ 2 — (s + a + b) x±, 2 +b = we obtain 



1_ f <+ M+1 (x7 m+L+1 - x 2 m+L+1 ) + 4 +M+1 (x 2 m+L+1 - V" +L+1 ) ,n<m 

nl r-n ~L+1 ( -m M+1 _ „-m„M+l\ , n „L+1 ( -m M+1 _ „-m„M+l\ > 
/ s 1 x t L ^ V 1 2 2 1 7 2 1 V 2 1 1 2 7 ) — 



where D = a {x x - x 2 ) {x{ +M+2 - x% +M+2 ) . Defining 1 f z - > 112 



the variable z = s + a + b — 2/j,, the constant fx = \/ah, B(z) ^ + ^ 2 j . (B5b) 
and the functions 

/ 2 \i/ 2 _ 

A(z) = l-\ - + — g , (B5a) results in x 2 = ^ A(z) and a(x 1 -x 2 ) = B(z) \ Using 

2fi \[i 4fj, J these functions we can recast the expression (B4), for the 
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case M = L, as 

G S nm (s) = l^B(z) [l-A(z) 4 < £+1 > 
x A(z)l™" m l 

-A{z) 2(L+1) ( K A{z) n+m + A(z)-( n+m 
-M(z) 4(L+1) A(z)^ n - m 



(a,b),(n,m) <-> (6, a),(m,n). Denoting by Gn m (s) the 
Green's function corresponding to the survival problem 
without bias (and taking /i = y/ab) we can im- 

mediately see that G^ m (s) = ^'^V 2 - Q^ m (z). Moreover, 
give n th at (z — /x/C°) G (z) = 1, we immediately obtain 
Eq. (IB!]). 



(B6) 



The last expression satisfies the symmetry relation 



We now compute the propagator J„ m (s) = K° G n m (s). 
Applying the operator K° defined in Eq. (3.3) to 
Eq. (pq), the following expression is found: 



J 



J% m (s) = B(z) [l-A(z) 4 ^ 



A{z) 



\^A{z)] r 



A(z) 



A(z) 2(L + l) 

V7 , A(z) 



A(z) 
VlA(z) 



V7 

1 



^A(z) 



where n,m € D and 

A{z) 

Z nm = { VlA{z 



V7 

-2 + A(z) 



(^A(z))-™- 1 (/7A(z)-l) 2 



A{z) yfyA(z) 



^ (^Aiz))— 1 {^A{z)-lf+A{zf^ (J^) 



V7 



- 1 



,n < m 



,n = m 



,n > m 



(B7) 



(B8) 



On the other hand, for the residence problem we pro- 
pose 



_ J Ai x™ for n < m 



G nm {s) <j A2 T „ for n > m 



with 



n, m efl. (B9) 



This expressi on immediately satisfies the homoge neou s 
case of Eq. (Bl) and the boundary conditions (B2t), 
given that x\ xi = 7 < 1 with x\ > 1. The constants 
Ai and A2 are fixed by imposing for n = m, the equal- 
ity b etw een both expressions in Eq. ( |B9| ), and writting 
Eq. (Bl) for n = m. Hence, we obtain 



1 J 7" ^2 
a(xi - x 2 ) L £ 2 



_ m 
for n> m . 

(BIO) 



Using the functions defined by Eq. (B5), the last expres- 
sion can be recast as 



G« m ( S )= 7 C«- m >/ 2 B(z)A(z)\ 



(Bll) 



Applying the operator K° to Eq. ( Bll ) we find: 

J%n(s) = 7 (n - m)/2 S( 2 )A(z)l»- m l 

r 7V2 A(z)- 1 + 7- 1 / 2 A(z) - 2 ,n < m 
x <| 7V2 4(3) + 7 -i/2 ^4( z ) _ 2 ,n = m 
{ 7V2 A(z) + 7-V2 A(z)- 1 - 2 ,n > m. 

(B12) 

We are interested in analyzing the propagator J nm (s) 
in the limit s — > for a given bias, also we are concerned 
with the expansion for small bias (e — > 0). We must stress 
that these limits do not commute. Thus, we first must 
take the limit s -* for a fixed bias, and then perform 
the expansion in the parameter e. Defining the variable 
y = s/(a — b), we get from Eq. (B5), for small s, that 



■ y 



(B13a) 



B(z) 



1 



o(l-7) 



1 _i±l, + (l+7) 2 + 27 / 



1-7 



(1-7) 2 



(B13b) 
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These approximations allow us to write the propagator 
in the form: J nm (s) w A nm + 0(y), where for the survival 
problem 



^ a (1-7) V 1 J 

(1 - 7 2 ) (-7 ,l+i + 7 ™-™+2i+i^\ 



, n < m 

> < -(1 - 7 2 j 7 "+ L + (l - 7 ) ( 7 2L +! - 1) , n = m 
(1 - 7 2 ) (- 7 "+ L + 7 «-™-i) , n > m, 

(B14) 

and, on the other hand, for the residence problem 




A R = — 

nm a 



1 



(1-7)7-' 



, n < m 
,n = rn 
,n > m 



(B15) 



We remark that we keep the exact dependence in the 
parameter L (s ize of the system) for the FPT problem. 
From Eq. ( B14 ), we can also see that the s-independent 
contribution is quite different from that which is obtained 
without bias |1(| . In the presence of bias we get nondiag- 
onal contributions to the propagator for the survival and 
residence problems in the limit s-*0, which also remain 
for small bias, as has been shown in Eq. (5.17). 



APPENDIX C: TERWIEL'S CUMULANTS 

The Terwiel's cumulants were introduced in Ref. p5| . 
In this Appendix we are concerned with the evaluation 
of the cumulants of the random operator tp n defined by 



where 



ijj n = M n {\-N n VM n ) 



M„ 



1 Jn 



is a random variable, and 

N n = 
From its definition, 



J mi 



1 - (M n ) J n 



(CI) 



(C2) 



(C3) 



{i>ni>mi ■ ■ ■ 4>m p )T = V t/j n (1 - V) Ipmi ■ ■ ■ (1 ~ V) 1pm p 

(C4) 

can be easily obtained 

{i>n)T = (ipn) , 

-(V'nV'miXV'ma) + ( < 0n)(V , rai)(V , m 2 ) • 

(C5) 



For n 7^ rn, ip n and i/j m are statistically independent, 
then (ip n ip m ) = (ip n )(4 > m)- Therefore, 

(i>ni>m)T = 0, 
(ipnlpm^n')T = {n ^ Tl') , 

(lpn4> m 1pn}T = (ipnlpmlpn) + (^n) (ipm) (i>n) > (C6) 



and using Eqs. (Cl)-(C3h results 



_ (M n ) 
{Pn) l + (M n )J nn > 

(tpntp m 1pn}T = (1 + (M n ) J„ n ) ^ 

X ((M 2 )(M m ) - N m (M n ) 2 (M m ) 2 ) 



(C7) 



In the limit s — > 0, from Eq. ( ^5.17 ) the diagonal com- 
ponents of the propagator J nn can be written as 



where 



--(i + M > 



^+Tiy , FPT problem 
, RT problem . 



Using the transfer rate oj = a + , we obtain 



(C8) 



(C9) 



(CIO) 



(CH) 



For a disorder of class (]a|) for which the quantities ft 
are finite, results in 



( 1 £,n Jnn) 

Thus, we can write 



1- (1-- 













l - 


('-=) 













(l-i)> = l-afl, 
(l--M = l-2aft+a 2 ft. (C12) 



Hence, 

(M„) f 
1 + (M„)J„„ : 

and therefore 



a[l-aft - (l-2aft +a 2 ft)0„e] , 



oft 
1 



1 - 1 -a- 



■ft" 



ft 



2 - a 



ft 



(C13) 



ft 1 — a — a J 7 8 n e , 

ft" 2 [a (1 - a ft ) (l - 2a ft + a 2 ft) 



(1- aft) 4 +0(e), 



(C14) 
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for n 7^ to. Here, we have used that T = (fc/Pf) — 1. Finally, for a disorder of class ([]), ft cx e 1 - fe - Q and re- 

For the strong classes of disorder, the quantities ft suits in (M n ) cx e~ a and (M^) cx e _1_ ". Therefore, 
diverge in the limit e — > and can be seen that 

1 V\ cx ft. In particular, for disorder of class (Q) 



-k 



(3% cx | lne| and ft a 6 1 1 for (k > 1). Thus, from 
Eq. ( pTT| ) we obtain (M„) cx |lne| and (M^) cx e _1 , oc ft -1 cx e a , 

which give (iMwMr oc e- 1 . (C16) 

<V^> ~ /Sf 1 oc ' 



In 



(lpnlf> m *Pn)T OC 1 . (C15) 

e me 
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